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Abstract 



We introduce the differential, integral, and variational delta-embeddings. We prove that the in- 
tegral delta-embedding of the Euler-Lagrange equations and the variational delta-embedding 
^^ , coincide on an arbitrary time scale. In particular, a new coherent embedding for the discrete 

^^ ■ calculus of variations that is compatible with the least action principle is obtained. 
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'nI" ' 1. Introduction 

Cn . An ordinary differential equation is usually given in differential form, i.e.. 



g I ^^f(t' -^(0) , te[a,b], x(t) e W. 

However, one can also consider the integral form of the equation: 

x(t) = x(a) -\- I f(s,x(s))ds, te[a,b]. 

Ja 



X 

d • The differential form is related to dynamics via the time derivative. The integral form is useful 

for proving existence and unicity of solutions or to study analytical properties of solutions. 

In order to give a meaning to a differential equation over a new set (e.g., stochastic processes, 
non-differentiable functions or discrete sets) one can use the differential or the integral form. In 
general, these two generalizations do not give the same object. In the differential case, we need 
to extend first the time derivative. As an example, we can look to Schwartz's distributions 1 15| or 
backward/forward finite differences in the discrete case. Using the new derivative one can then 
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generalize differential operators and then differential equations of arbitrary order In the integral 
case, one need to give a meaning to the integral over the new set. This strategy is for example used 



by Ito 111 in in order to define stochastic differential equations, defining first stochastic integrals. 
In general, the integral form imposes less constraints on the underlying objects. This is already 
in the classical case, where we need a differentiable function to write the differential form but 
only continuity or weaker regularity to give a meaning to the integral form. 

The notion of embedding introduced in |7] is an algebraic procedure providing an extension 
of classical differential equations over an arbitrary vector space. Embedding is based on the 
differential formulation of the equation. This formalism was developed in the framework of 
fractional equations |6], stochastic processes |7], and non-differentiable functions |8]. Recently, 
it has been extended to discrete sets in order to discuss discretization of differential equations 
and numerical schemes ||5t] . 

In this paper, we define an embedding containing the discrete as well as the continuous case 
using the time scale calculus. We use the proposed embedding in order to define time scale 
extensions of ordinary differential equations both in differential and integral forms. 

Of particular importance for many applications in physics and mathematics is the case of 
Lagrangian systems governed by an Euler-Lagrange equation. Lagrangian systems possess a 
variational structure, i.e., their solutions correspond to critical points of a functional and this 
characterization does not depend on the coordinates system. This induces strong constraints on 
solutions, for example the conservation of energy for autonomous classical Lagrangian systems. 
That is, if the Lagrangian does not depend explicitly on t, then the energy is constant along 
physical trajectories. We use the time scale embedding in order to provide an analogous of the 
classical Lagrangian functional on an arbitrary time scale. By developing the corresponding cal- 
culus of variations, one then obtain the Euler-Lagrange equation. This extension of the original 
Euler-Lagrange equation passing by the time scales embedding of the functional is called a time 
scale variational embedding. These extensions are known under the terminology of variational 
integrators in the discrete setting. 

We then have three ways to extend a given ordinary differential equation: differential, inte- 
gral or variational embedding. All these extensions are a priori different. The coherence problem 
introduced in |7[|, in the context of the stochastic embedding, consider the problem of finding 
conditions under which these extensions coincide. Here we prove that the integral and varia- 
tional embeddings are coherent (see Theorem [TT}. The result is new and interesting even in the 
discrete setting, providing a new form of the Euler-Lagrange difference equation (see (fTSl l) that 
is compatible with the least action principle. 

2. Note on the notation used 

We denote by / or f — > f{t) a function, and by f{t) the value of the function at point t. Along 
all the text we consistently use square brackets for the arguments of operators and round brackets 
for the arguments of all the other type of functions. Functionals are denoted by uppercase letters 
in calligraphic mode. We denote by D the usual differential operator and by 5, the operator of 
partial differentiation with respect to the ith variable. 

3. Reminder about the time scale calculus 

The reader interested on the calculus on time scales is refereed to the book f4\. Here we just 
recall the necessary concepts and fix some notations. 
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A nonempty closed subset of M. is called a time scale and is denoted by T. Thus, M, Z, 
and N, are trivial examples of time scales. Other examples of time scales are: [-2,4] U^^^ 
hZ :- {hz\z e Z) for some h > 0, q^" :- {q''\k e No] for some ^ > 1, and the Cantor set. 
We assume that a time scale T has the topology that it inherits from the real numbers with the 
standard topology. 

The forward jump cr : T — > T is defined by cr{t) = inf {i 6 T : 5 > f} for all f e T, while 
the backward jump p : T — > T is defined by p{t) = snp{s eT : s < t} for all f e T, where 
inf = sup T (i.e., cr{M) = M if T has a maximum M) and sup = inf T (i.e., pirn) = m if T has 
a minimum m). The graininess function /i : T — > [0, oo) is defined by ^{t) - cr{t) - 1 for all f E T. 

Example 1. IfT - E, then cr{t) - p{t) - tandid(t) - 0. IfT - hZ, then cr{t) -t + h, p(t) - t-h, 
and fj.(t) — h. On the other hand, ifT — q^°, where q > \ is a fixed real number, then we have 
cr(t) - qt, p(t) - q'^t, and fi{t) - iq - l)t. 

In order to introduce the definition of delta derivative, we define a new set T* which is derived 
from T as follows: if T has a left-scattered maximal point M, then r := T \ {M}; otherwise, 
T* := T. In general, for r > 2, T* \- ST" \ . Similarly, if T has a right-scattered minimum m, 
then we define T,^ := T \ \m\; otherwise, T^^ \- T. Moreover, we define T* := T n T^^. 

Definition 1. We say that a function / : T — > K /i delta differentiable att^f^ if there exists a 
number A[/](f) such that for all e > there is a neighborhood U oft such that 

\f(cr(t)) - f(s) - A[f](t)(cr(t) - s)\ < s\cr(t) - s\, for all s e U. 

We call A[/](f) the delta derivative of f at t and we say that f is delta differentiable on T" 
provided A[/](f) exists for all t e T". 

Example 2. IfT — M, then A[f](t) — D[f](t), i.e., the delta derivative coincides with the usual 
one. IfT - liL, then A[/](f) = j{f{t + h) — f{t)) —: A+[/](f), where A+ is the usual forward 
difference operator defined by the last equation. IfT — ^^°, q> \, then A[/](f) = i _()] , i-^-, 
we get the usual derivative of quantum calculus. 

A function / : T — > M is called rd-continuous if it is continuous at right-dense points and 
if its left-sided limit exists at left-dense points. We denote the set of all rd-continuous functions 
by Crd(T, K) and the set of all delta differentiable functions with rd-continuous derivative by 
Cj'j(T, M). It is known (see ||4|, Theorem 1.74]) that rd-continuous functions possess a delta 
antiderivative, i.e., there exists a function f with A[^] - /, and in this case the delta integral is 
defined by // f(t)At = ^(d) - ^(c) for all c,deT. 



Example 3. Let a,b eT with a < b. IfT = M, then J fit) At - J f(t)dt, where the integral on 

rh --1 

the right-hand side is the classical Riemann integral. IfT — /iZ, then J f(t)At — Yi'l-i hf(kh). 
//T = ^N„^ ^ > 1^ f^^^ Jb ^^^^^^ ^(l-q) 2,^j^ ^, tf(t). 

The delta integral has the following properties: 

(i) if / e Cd and t e T, then 



I 



■o-(l) 

/(t)At = M0/(0; 



(ii) if c,d e T and / and g are delta-differentiable, then the following formulas of integration 
by parts hold: 

J f(cr(t))MgmAt = (fgmfrj - j M/mgiOAt, 
J mMgmAt = (fgmfrJ - J Ai/mgicrimAt. 



(1) 



4. Time scale embeddings and evaluation operators 

Let T be a bounded time scale with a := min T and b :- max T. We denote by C([a, b];M.) 
the set of continuous functions x : [a,b] — > M. As introduced in Section |3] by CrdCT, M) we 
denote the set of all real valued rd-continuous functions defined on T, and by C^.^(T, M.) the set 
of all delta differentiable functions with rd-continuous derivative. 

A time scale embedding is given by specifying: 

• A mapping i : C([a, b],W) ^ C,-d(T, R); 

• An operator 6 : C^{[a, b],W) — > CJ.^{V, K), called a generalized derivative; 

• An operator J : C{[a, b], M) —> Crd(T, M), called a generalized integral operator. 

We fix the following embedding: 

Definition 2 (Time scale embedding). The mapping l is obtained by restriction of functions to 
T. The operator 6 is chosen to be the A derivative, and the operator J is given by the A integral 

as follows: 

ro-(i) 

6[x](t):^A[x](t), J[x](t):^ x(s)As . 



Definition 3 (Evaluation operator). Let / : R — > M fee a continuous function. We denote by f the 
operator associated to f and defined by 



7: ^(^'^^ -- ~ — (2) 

^ X ^ f[x] := t ^ f(x(t)). ^ ' 

The operator f given by (|2]l is called the evaluation operator associated with f. 

The definition of evaluation operator is easily extended in various ways. We give in Defini- 
tion |4] a special evaluation operator that naturally arises in the study of problems of the calculus 
of variations and respective Euler-Lagrange equations (cf. Section|5]|. 

Definition 4 (Lagrangian operator). Let L : [a, fe] x K x M — > M fee a C' function defined for all 
it, X, v) 6 [a, fe] X R2 by Lit, x, v) e R. The Lagrangian operator L : 0(\a, fe], R) -^ C'([a, b\W) 
associated with L is the evaluation operator defined by L\x\ \— t ^> L(t, jc(f), D\x\(t)). 



We consider ordinary differential equations of the form 

O[x](f) = 0, te[a,b], 
where x e C"(M., M) and O is a differential operator of order n,n> 1, given by 



= Yj^r{D'ob), (3) 

(=0 

where (ai) (resp. (/?,)) is the family of evaluation operators associated to a family of functions 
(fl/) (resp. (fe,)), and D' is the derivative of order /, i.e., D' - jp. Differential operators of form 
([3]l play a crucial role when deahng with Euler-Lagrange equations. 

We are now ready to define the time scale embedding of evaluation and differential operators. 

Definition 5 (Time scale embedding of evaluation operators). Let f : M. -^ M. be a continuous 
function and f the associated evaluation operator. The time scale embedding fj of f is the 
extension off to C,,i(T, E).- 

~ C-rf(T,M) -^ __ C,,;(T,M) 
■"■ X ^ fT[x] -.^ t ^ f(x(t)). 

Next definition gives the time scale embedding of the differential operator (l3)- 

Definition 6 (Time scale embedding of differential operators). The time scale embedding of the 
differential operator (O is defined by 

n 

Oa = ^ aij ■ (A' o bij) . 

1=0 

The two previous definitions are sufficient to define the time scale embedding of a given 
ordinary differential equation. 

Definition 7 (Time scale embedding of differential equations). The delta-differential embedding 
of an ordinary differential equation 0[x] — 0, x e C"{[a,b],W), is given by Oa[x] — 0, x & 

C"JT",R). 

In order to define the delta-integral and variational embeddings (see Sections |2l[l] and |9]l we 
need to know how to embed an integral functional. 

Definition 8 (Time scale embedding of integral functionals). Let L : [a,b] xM?' -^ M. be a 
continuous function and £, the functional defined by 

Z{x)^ I L{s,x{s),D[x](s))ds^ I L[x]{s)ds. 

The time scale embedding -£a of £. is given by 

£a(x)- I L(s,x(s),A[x](s))As- I Lt[x](s)As . 
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5. Calculus of variations 

The classical variational functional £. is defined by 

w= r 

J a 



■b 

£(x)^ I L(t,x(t),D[x](t))dt, (4) 



where L : [a, &] X R x M — > K is a smooth real valued function called the Lagrangian (see, 
e.g., ri61)- Functional (|4| can be written, using the Lagrangian operator L (Definition |4|, in the 
following equivalent form: 

£(x) = I L[x](t)dt. 
The Euler-Lagrange equation associated to (|4) is given (see, e.g., Il16J) by 

D[t^ di[L] (t, x(t), D[x](t))] it) - d2[L] it, x{t\ D[x]{t)) = 0, (5) 

t e{a,b'\, which we can write, equivalently, as 

[d o d^L]) [xm - dAL][xm = 0. 

Still another way to write the Euler-Lagrange equation consists in introducing the differential 
operator EL/,, called the Euler-Lagrange operator, given by 

ELi := D o d^] - dAL] . 
We can then write the Euler-Lagrange equation simply as EL£[;t:](f) = 0, f e [a, b]. 

6. Delta-differential embedding of the Euler-Lagrange equation 

By Definition |6] the time scale delta embedding of the Euler-Lagrange operator EL/, gives 
the new operator 

(ELJa := A o {dAL])^ - {WX ■ 
As a consequence, we have the following lemma. 

Lemma 9 (Delta-differential embedding of the Euler-Lagrange equation). The delta-differential 
embedding of the Euler-Lagrange equation is given by (ELi)iAx\{t) = 0, f e V , i.e.. 



(a o drnj) [xm - drnjixm = o (6) 



for any t eV 



In the discrete case T - {a, b] n hZ, we obtain from (|6]l the well-known discrete version of 
the Euler-Lagrange equation, often written as 

dL dL 

A+ o — (f, x(t), A+x(0) - -TT (t, xit), A+x(0) = 0, (7) 

ov ox 

t € T" , where A+/(f) = liilipM xhe important point to note here is that from the numerical 
point of view, equation Q does not provide a good scheme. Let us see a simple example. 
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Example 4. Consider the Lagrangian Lit, x, v) = j v^ - U{x), where U is the potential energy 
and (t, X, v) e [a, fe] X R X R. Then the Euler— Lagrange equation (|7]l gives 

Xk+2 - 2xk+i + Xk dU 

72 '^'dl^^''^" k^0,...,N-2, (8) 

where N = ^ and x^ = x{a + hk). This numerical scheme is of order one, meaning that we 
make an error of order h at each step, which is of course not good. 

In the next section we show an ahemative Euler-Lagrange to (|7]) that leads to more suitable 
numerical schemes. As we shall see in Section |9] this comes from the fact that the embedded 
Euler-Lagrange equation (|6) is not coherent, meaning that it does not preserve the variational 
structure. As a consequence, the numerical scheme (jH) is not symplectic in contrast to the flow 
of the Lagrangian system (see 1 1]). In particular, the numerical scheme ([8]) dissipates artificially 
energy (see [.13. Fig. 1, p. 364]). 

7. Discrete variational embedding 

The time scale embedding can be also used to define a delta analogue of the variational 
functional (|4|i. Using Definition[8] and remembering that cr{b) = b, the time scale embedding of 
dill is 

XaW= L(f,x(0,A[x](f)) Af= Lf[x](t)ht. (9) 

A calculus of variations on time scales for functionals of type ^ is developed in Section|9] Here 
we just emphasize that in the discrete case T - [a, b~\ n /iZ functional (|9]l reduces to the classical 
discrete Lagrangian functional 

N-\ 

£A(x)^hY^L(tk,Xk,A+Xk), (10) 

where A^ = ^, x^ - x(a+hk) and A+x^ - "''^^^~"''' , and that the Euler-Lagrange equation obtained 
by applying the discrete variational principle to (ITOt takes the form 

A_ o -^ (t, x(t), A+x(0) - -^ (t, x(t), A+x(0) = 0, (11) 

ov ox 

t € TJ^, where A_ is the backward finite difference operator defined by A_/(f) = ^ f5,'l0tl- 
The numerical scheme corresponding to the discrete variational embedding, i.e., to ( fTTT ). is 
called in the literature a variational integrator IsL llOll . Next example shows that the variational 
integrator associated with problem of Example 2]is a better numerical scheme than dS). 

Example 5. Consider the same Lagrangian as in Example^of Section^ L(t, x,v) — \y^ — U{x), 
where (t, x, v) e [a,b] xWx M. The Euler-Lagrange equation (llll i can be written as 

Xk+\ - 2xk + Xk-\ dU 

1? "^57^^*-^" k^l,...,N-\, 

where N — -jf- and Xk = x(a + hk). This numerical scheme possess very good properties. In 
particular, it is easily seen that the order of approximation is now two and not of order one as in 
Example^ 

7 



We remark that the form of Xa given by ^ is not the usual one in the hterature of time scales 
(see II2L I12II14I1 and references therein). Indeed, in the literature of the calculus of variations on 
time scales, the following version of the Lagrangian functional is studied: 



Ja 



j:^^"-^\x)= L(t,x(o-(t)),A[x](t))At. (12) 

Ja 

However, the composition of x with the forward jump cr found in (fT2l i seems not natural from 
the point of view of embedding. 

8. Delta-integral embedding of the Euler-Lagrange equation in integral form 

We begin by rewriting the classical Euler-Lagrange equation into integral form. Integrating 
(|5]i we obtain that 



Ja 



d3[L](t,x{t),D[xm)= d2[L](T,x(T),D[x](T))dT + c, (13) 

Ja 

for some constant c and all t e [a,b] or, using the evaluation operator, 

dAL][xm = r dAL][x]{T)dT + c. 



Using Definition[8] we obtain the delta-integral embedding of the classical Euler-Lagrange equa- 
tion ( fBT l. 

Lemma 10 (Delta-integral embedding of the Euler-Lagrange equation in integral form). The 
delta-integral embedding of the Euler-Lagrange equation (I13l l is given by 

d3[L](t,x(t),A[x]it))^ d2[L](T,x(T),A[x](T))AT + c (14) 

Ja 

or, equivalently, as 

']tw(o= r 

Ja 



53[L]tW(0= I d2[Lh[x]{T)AT + c. 



where c is a constant and t e T*^. 

Note that in the particular case T = {a, b} Pi /iZ equation (fT4l i gives the discrete Euler- 
Lagrange equation 

dL I Xk+\ - Xk\ -fT^ dL I Xi+i - X, 



I Xk+i-Xk\ ^ dL I Xi+i-Xi\ 

(..,x,,^— ) = /.^-(,,x„^— ) + c, (15) 



cJv . ... 

(-0 



where ti — a + hi, i - 0, . . .,k, x/ - x(f,), and A: = 0, . . . , A^ - 1 . This numerical scheme is different 
from O and ( fTTT i. and has not been discussed before in the literature with respect to embedding 
and coherence. This is done in Section |9] 



9. The delta-variational embedding and coherence 

Our next theorem shows that equation (fT4l i can be also obtained from the least action prin- 
ciple. In other words, Theorem [TT] asserts that the delta-integral embedding of the classical 
Euler-Lagrange equation in integral form (ITji and the delta-variational embedding are coherent. 

Theorem 11. If x is a local minimizer or maximizer to ^ subject to the boundary conditions 
x{a) — Xa and x(b) — Xh, then x satisfies the Euler-Lagrange equation ( 1141 ) for some constant c 
and all t e T". 



Proof. Suppose that Xa has a weak local extremum at x. Let x - x + eh, where e e M is a small 



parameter, h e Cy, such that h{a) = h{b) - 0. We consider 



0(e) := £a{x + sh) =1 L(t, x(t) + sh(t), A[x](t) + sA[h](t)) At. 

J a 

A necessary condition for x to be an extremizer is given by 

0'(e)|^^„ = o {d2[L] it,x(t),A[x](t))h(t) + d^L] (t,x(t),A[xm)A[hm)At = 0. (16) 

Ja 

The integration by parts formula ([1) gives 



/ 



■b 

d2[L](t,m,A[x]{t))Kt)At 



52[L](t,x(t),A[x](t))At/z(0I;:^ - J J d2[mT,x{T),A{x\(r))ATA[h]{t)\At. 



Because h{a) - h(b) - 0, the necessary condition (fTSI l can be written as 

Xb j r-a-(t) \ 

\dAL](t,x(t),A[x\{t))- I d2[L\{t,x{t\A[x\{t))AT\A[h]{t)At = 

for all h 6 C].^^ such that h{a) - h(b) - 0. Thus, by the Dubois-Reymond Lemma (see yi 
Lemma 4.1]), we have 

dilL] it, m, A[x](0) = d2[L] (t, xiT), A[Jc](t)) At + c 

Ja 

for some c e M and all teT". D 



10. Conclusion 

Given a variational functional and a corresponding Euler-Lagrange equation, the problem of 
coherence concerns the coincidence of a direct embedding of the given Euler-Lagrange equa- 
tion with the one obtained from the application of the embedding to the variational functional 
followed by application of the least action principle. An embedding is not always coherent and 
a nontrivial problem is to find conditions under which the embedding can be made coherent. 
An example of this is given by the standard discrete embedding; the discrete embedding of the 
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Euler-Lagrange equation gives (|7} but the Euler-Lagrange equation (fTTl l obtained by the stan- 
dard discrete calculus of variations does not coincide. On the other hand, from the point of 
view of numerical integration of ordinary differential equations, we know that the discrete vari- 
ational embedding is better than the direct discrete embedding of the Euler-Lagrange equation 
(cf. Example |5]l. The lack of coherence means that a pure algebraic discretization of the Euler- 
Lagrange equation is not good in general, because we miss some important dynamical properties 
of the equation which are encoded in the Lagrangian functional. A method to solve this default of 
coherence had been recently proposed in [5J, and consists to rewrite the classical Euler-Lagrange 
(|5]l as an asymmetric differential equation using left and right derivatives. Inspired by the results 
of ||9[], here we propose a completely different point of view to embedding based on the Euler- 
Lagrange equation in integral form. For that we introduce the new delta-integral embedding (cf. 
Definition [8]l. Our main result shows that the delta-integral embedding and the delta- variational 
embedding are coherent for any possible discretization (Theorem[TT]is valid on an arbitrary time 
scale). 
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